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Abstract.  Direct  numerical  simulation  of  coagulating  aerosols  in  two- 
dimensional,  mixing  layers  are  performed.  The  evolution  of  the  particle 
field  is  obtained  by  utilizing  a sectional  model  to  approximate  the  aerosol 
general  dynamic  equation.  The  sectional  model  is  advantageous  in  that 
there  are  no  a priori  assumptions  regarding  the  particle  size  distribution. 
This  representation  facilitates  the  capture  of  the  underlying  physics  in  an 
accurate  manner. 


1.  Introduction 

Ultrafine  particles  play  a very  important  role  in  a wide  variety  of  phys- 
ical/chemical phenomena  and  processes.  An  important  application  is  the 
synthesis  of  nanostructured  materials.  There  are  several  technologies  which 
can  be  employed  in  the  manufacture  of  nanoscale  materials  (films,  parti- 
cles, etc.)  Vapor-phase  methodologies  are  by  far  the  most  favored  because 
of  chemical  purity  and  cost  considerations.  A key  issue  in  the  formation 
of  nanoscale  particles  is  the  prevention  of  hard  agglomerates  and  chemical 
control.  A number  of  strategies  have  been  attempted  to  minimize  agglom- 
eration(Matsoukas  and  Friedlander,  1991). 
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The  dynamics  of  particles  in  turbulent  flows  have  received  some  atten- 
tion. The  earlier  works  were  focused  on  understanding  the  phenomenon  of 
particle  dispersion  by  turbulence(Riley  and  Patterson,  1974;  Elghobashi, 
1991).  The  influence  of  particle  parameters  on  collision  frequencies  in  a 
turbulent  particle  laden  suspension  leading  to  coagulation  was  considered 
by  Sundaram  and  Collins  (1996),  who  showed  that  the  magnitude  of  the 
minimum  particle  collision  frequency  was  strongly  correlated  with  the  tur- 
bulent motions  at  the  integral  scale.  Reade  and  Collins  (2000)  simulated 
the  coagulation  and  growth  of  an  initially  mono-disperse  aerosol  subject  to 
isotropic  turbulence.  This  work  resulted  in  an  improved  understanding  of 
the  trends  in  the  relative  width  of  the  particle  size  distribution  and  its  de- 
pendence on  the  Stokes  number  and  radial  distribution  function.  However, 
much  of  the  work  performed  thus  far  consider  large,  micron-scale  particles 
in  a Lagrangian  manner  which  use  primarily  particle  tracking  methods.  The 
large  number  of  particles  needed  to  represent  the  underlying  physics  of  par- 
ticle growth  using  Lagrangian  methods  render  the  computations  infeasible 
for  all  but  inhomogeneous  systems. 

In  this  work  we  consider  the  growth  of  nanoscale  particles  in  a com- 
pressible, temporally  developing  mixing  layer.  The  particulate  phase  will 
be  accounted  for  using  a sectional  method  which  treats  the  particles  in 
an  Eulerian  manner (Pyykonen  and  Jokiniemi,  1999;  Garrick  et  al.,  2001). 
This  approach  is  advantageous  in  that  there  are  no  a priori  assumptions  re- 
garding the  nature  of  the  particle  size  distribution.  We  intend  to  use  direct 
numerical  simulations  to  capture  the  underlying  physics  of  particle  growth 
in  a model- free  manner  (Givi,  1989). 

2.  Formulation 

The  flows  under  consideration  are  two-dimensional  mixing  layers  and  are 
governed  by  the  compressible  Navier-Stokes  equations.  The  transport  of  the 
nanoscale  particles  dispersed  is  governed  by  the  aerosol  general  dynamic 
equation  (GDE).  The  equation  is  expressed  in  a discrete  form  as  a popu- 
lation balance  on  each  cluster  or  particle  size.  Prom  a practical  standpoint 
however  such  systems  of  equations  cannot  be  solved  explicitly  except  for 
very  small  particle  sizes.  Therefore  a sectional  method  is  used  to  represent 
the  particle  field  as  a function  of  space  and  time. 

This  approach  effectively  divides  the  particle  size  distribution  into  “bins,” 
as  illustrated  in  Fig.  1.  The  discrete  part  of  the  representation  models  nucle- 
ation  and  molecule-molecule  interactions  which  lead  to  particle  formation. 
A molecular  cluster  qi  is  comprised  of  i molecules.  Typically  molecular 
clusters  grow  by  the  addition  of  one  molecular  unit  at  a time.  As  parti- 
cles become  larger  they  are  transitioned  to  the  sectional  representation, 
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Figure  1.  Discrete  - sectional  representation  of  particle  field. 


Qk-,  after  reaching  a cut-off  size  consisting  of  M clusters.  We  consider  only 
large  clusters,  and  particles  which  typically  contain  tens  of  thousands  of 
molecules  for  which  a discrete  representation  is  unnecessary.  The  GDE  is 
therefore  solved  as  a set  of  Ns  transport  equations,  one  for  each  section  Qk, 
k = 1,2,...,  Ns  (Xiong  and  Pratsinis,  1993).  In  adopting  this  framework  we 
can  write  the  general  transport  equation  for  the  concentration  of  particles 
in  the  Arth  section,  Qk- 

dpQk  dpQkUj 

dt  dxj 

where  Vq  is  the  diffusivity 


^“^37 r/<  w 

where  kf,  is  the  Boltzmann  constant,  Cc  is  the  Cunningham  correction  fac- 
tor, dp  is  the  particle  diameter,  and  T and  fi  are  the  fluid  temperature  and 
viscosity,  respectively.  The  source  term,  iv®,  is  given  by 


9Qk 

dxj 


+ cp 


Q 
k > 


(1) 


Cc 


1 Ns  Ns  Ns 

“f-  EE  Pi,j  XijkQiQj  'Y^PikQiQk-  (3) 

i=  1 j= 1 i=  1 

The  source  term,  represents  the  effects  of  particle-particle  interactions: 
production  of  Qk  due  to  collisions  of  smaller  particles;  the  loss  or  gain  of  Qk 
by  collision  with  a smaller  particle  which  either  moves  the  resulting  particle 
out  of  or  into  section  k ; the  loss  of  particles  in  section  k as  they  collide  with 
each  other  and  form  larger  particles;  and  the  loss  of  particles  in  section  k due 
to  collisions  with  larger  particles.  It  should  be  noted  that  repeated  indices 
in  Eq.  3 do  not  imply  summation  but  instead  infer  interactions  between 
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particles  in  section  i and  particles  in  section  j.  The  collision  frequency 
function  is  given  by 


(4) 


where  Vi  is  the  volume  of  the  ith  particle,  pp  is  the  particle  density  and  Xijk 
is  given  by 


Xijk  — N 


( Vk -H-pJj+y,-) 
Vk+l-Vk 
(Vi+Vj)-Vk-1 


Vk-Vk-l 


0 


if  Vk  < Vi  + Vj  < vk+i 
if  Vk- 1 < Vi  + Vj  < Vk 
otherwise. 


(5) 


The  sectional  method  is  discretized  in  size  space  such  that  the  volume  of 
particles  in  two  successive  sections  is  doubled,  i.e.  vk  = 2 x Vk-\.  This 
scheme  allows  us  to  span  a volume  range  of  V = vi  to  V = 2N“~l  x v\. 


3.  Numerical  Procedure 

The  governing  transport  equations  are  solved  using  a hybrid  MacCormack 
based  compact  difference  scheme(Carpenter,  1990;  Kennedy  and  Carpen- 
ter, 1994).  The  numerical  scheme  used  is  based  on  the  one-parameter  family 
of  dissipative  two-four  schemes  (DCPS)  (Carpenter,  1990).  The  accuracy  of 
the  scheme  is  second  order  in  time,  and  fourth  order  in  space.  All  calcula- 
tions are  performed  on  a rectangular  uniformly  spaced  grid.  The  computa- 
tional algorithm  uses  the  message  passing  interface  to  execute  on  parallel 
computing  platforms.  The  treatment  of  the  particulate  phase  in  a Eulerian 
manner,  as  opposed  to  a Lagrangian  one,  helps  to  keep  processor-processor 
communication  to  a minimum. 


4.  Results 

The  flows  under  consideration  are  two-dimensional,  compressible,  shear  lay- 
ers. A schematic  is  shown  in  Fig.  2.  Periodic  boundary  conditions  are 
used  in  the  streamwise  direction,  while  the  free  shear  condition  is  im- 
posed in  the  cross-stream  direction.  The  gas  is  air,  and  the  velocity  is 
initialized  with  a hyperbolic  tangent  in  the  cross-stream  direction.  The 
Reynolds  number  is  based  on  the  vorticity  thickness,  5W,  and  the  veloc- 
ity difference,  AU  = U\  — U2,  = -AT<^. . Two  cases  are  considered. 

Case  (I),  considers  particle  coagulation  in  a uniform  temperature  flow  at 
Ti  = T2  — 300K.  Case  (II),  considers  the  flow  where  the  temperature 
of  the  particle-free  stream  is  T\  = 300K,  and  the  temperature  of  the 
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Figure  2.  Flow  configuration 


particle-laden  stream  is  T2  = 900  A".  Thermophoresis  effects  are  not  in- 
cluded in  these  calculations,  but  inclusion  of  such  effects  pose  no  addi- 
tional mathematical  or  computational  difficulties(Pratsinis  and  Kim,  1989; 
Pyykonen  and  Jokiniemi,  1999).  Both  cases  consider  an  initial  volume  frac- 
tion of  Vf  — 9.4  x 10  8.  The  volume  fraction  is  defined  as  the  ratio  of  the 
volume  occupied  by  the  particulate  phase  to  that  occupied  by  the  fluid.  A 
total  of  ten  sections  are  solved,  i.e.  Ns  = 10.  This  allows  for  the  solution 
of  particles  covering  a range  of  3 orders  of  magnitude  in  volume. 

Computations  are  performed  on  a domain  of  27t  x 27t  in  the  streamwise, 
x,  and  cross-stream,  y,  directions  on  15002  grid  points.  With  this  resolu- 
tion, flows  with  a Reynolds  number  of  Regu  = 200  are  well  resolved.  All 
calculations  were  performed  on  a SGI-Origin3800  supercomputer,  and  each 
calculation  simulated  up  to  a nondimensional  time  of  t*  = 14.25. 

Particle  concentration  profiles  obtained  from  Case  (I)  are  shown  in  Fig. 
3.  Cross-stream  profiles  of  Q 1,  Q2,  Q 7,  and  Q 10,  are  obtained  by  averaging 
in  the  streamwise  x-direction,  at  four  times.  All  values  of  Q are  normalized 
by  the  initial  number  of  particles  in  the  particle-laden  stream,  Q\0.  At  time 
t*  = 0,  Qi/Qio  = 1 in  the  particle-laden  stream,  and  Q\/Q\0  = 0 in  the 
particle  free  stream.  Figure  3a  indicates  the  decrease  from  the  initial  Q\ 
concentration  in  the  particle-laden  stream,  and  increase  in  concentration  in 
the  initially  particle-free  stream.  As  the  particles  collide,  they  coagulate  to 
form  larger  particles,  thereby  moving  out  of  section  1.  The  increase  in  par- 
ticles in  the  y > 0 region  is  due  to  dispersion.  As  the  mixing  layer  evolves 
the  particle-laden  stream  is  mixed,  via  convection,  with  the  particle-free 
fluid.  The  figure  also  reveals  a peak  in  Q 1 near  the  interface  of  the  two 
streams,  the  magnitude  of  which  decreases  with  time.  The  peak  is  due  to 
the  gradient  in  Q\  in  the  shear  region  which  results  in  a lower  growth  rate, 
in  comparison  to  that  observed  in  the  freestream  of  the  initially  particle- 
laden stream.  A similar  trend  is  observed  in  section  2,  shown  in  Fig.  3b. 
Particles  collide  and  coagulate  to  form  larger  particles  which  can  no  longer 
be  accommodated  in  Q2.  The  decrease  of  Q 2 in  time  is  accompanied  by 
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Figure  3.  Cross  stream  variation  of  particle  concentrations,  Case  (I):  (a)  Q i;  (b)  Q>\ 
(c)  Qr,  (d)  Qio. 


increases  in  the  higher  numbered  sections,  Qk , k = 3, 4, . . . , 10.  The  con- 
centration of  particles  in  section  7 is  shown  in  Fig.  3c.  The  trend  is  different 
from  that  observed  in  the  first  two  sections  in  that  the  maximum  concen- 
tration is  observed  in  the  initially  particle-laden  stream.  The  concentration 
increases  to  Qj/Qig  = 0.0044,  at  time  t*  = 6.01,  then  begins  to  decrease. 
However,  the  concentration  in  the  initally  particle- free  stream  increases  over 
the  same  period  of  time.  Near  the  interface  of  the  two  streams  there  is  some 
overlap  of  the  profiles  at  later  times,  t*  = 6.01, 9.13,  and  12.3.  This  “equilib- 
rium” may  be  attributed  to  the  mixing  of  the  particle-free  stream  with  the 
particle-laden  stream  as  the  vortex  develops.  The  concentration  of  particles 
in  section  10  is  shown  in  Fig.  3d.  The  figure  reveals  that  the  concentration 
Qio  increases  at  each  y location  across  the  mixing  layer.  There  are  fewer 
Q io  particles  in  the  core  of  the  vortex  in  comparison  to  the  initially  particle- 
laden stream.  This  further  reflects  the  reduced  growth  rate  observed  in  the 
first  two  sections.  Additionally,  the  concentrations  of  all  sections  are  seen 
to  spread  out  into  the  particle-free  stream  with  time  as  the  flow  develops. 
A quantitative  view  of  the  particle  field  for  the  isothermal  flow,  Case  (I), 
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Figure  4 • Instantaneous  particle  concentrations  contours,  Case  (I):  (a)  Qi;(b)  Q 2;  (b) 
Qi\  (b)  Q10. 


is  presented  in  Fig.  4.  This  figure  shows  instantaneous  contours  of  the  con- 
centrations in  sections  1,  2,  7,  and  10,  taken  at  time  t*  = 12.3.  In  addition 
to  spatial  concentration  variation,  the  diffusion  effects  are  also  evident.  The 
striation  thickness  decreases  as  the  particle  size  increases.  This  is  because 
larger  particles  have  smaller  coefficients  of  diffusion. 


Figure  5.  Cross  stream  variation  of  mean  particle  diameter:  (a)  Case  (I);  (b)  Case  (II). 
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Temperature  effects  may  be  observed  by  considering  the  average  particle 
size.  The  average  diameter  is  given  by  dp  = (6/7 rw)1/3,  where  the  mean 
volume  is  given  by  v = Qivi/Yli^  1 Qi-  The  temporal  evolution  of  the 

mean  diameter  for  Case  (I),  and  Case  (II)  is  shown  in  Fig.  5.  Though 
the  profiles  show  an  overall  increase  with  time,  the  particle  growth  rate 
in  Case  (II)  is  less  than  that  in  Case  (I).  Though  the  temperature  of  the 
particle-laden  stream  is  greater,  which  also  means  that  the  rate  of  collisions 
is  greater,  the  particles  effectively  grow  at  a slower  rate.  The  spread  of  the 
profiles  indicates  mixing  of  the  two  streams  as  the  vortices  develop. 
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